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RANDOM AND FREE POSITIVE MAPS WITH APPLICATIONS TO 

ENTANGLEMENT DETECTION 

BENOiT COLLINS, PATRICK HAYDEN, AND ION NECHITA 


Abstract. We apply random matrix and free probability techniques to the study of linear maps 
of interest in quantum information theory. Random quantum channels have already been widely 
investigated with spectacular success. Here, we are interested in more general maps, asking only for 
fc-positivity instead of the complete positivity required of quantum channels. Unlike the theory of 
completely positive maps, the theory of fc-positive maps is far from being completely understood, 
and our techniques give many new parametrized families of such maps. 

We also establish a conceptual link with free probability theory, and show that our constructions 
can be obtained to some extent without random techniques in the setup of free products of von 
Neumann algebras. 

Finally, we study the properties of our examples and show that for some parameters, they are 
indecomposable. In particular, they can be used to detect the presence of entanglement missed by 
the partial transposition test, that is, PPT entanglement. As an application, we considerably refine 
our understanding of PPT states in the case where one of the spaces is large whereas the other one 
remains small. 
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1. Introduction 

Completely positive maps play a privileged role in quantum information theory because they 
provide the quantum mechanical generalization of the notion of a noisy channel. The Stinespring 
theorem provides a useful classification of these maps that also justifies their physical relevance 
as ignorance-induced noisy summaries of the unitary evolution on a larger system [30]. Random 
methods have played an important role in the study of noisy quantum channels, notably in the proofs 
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of capacity theorems [27, 21, 17] and in the construction of channels with unusual properties [18, 
28, 14], including channels violating the additivity conjecture [19, 16, 12, 6]. 

Maps with weaker positivity properties are of interest in quantum information theory precisely 
because they fail in subtle ways to be physically realizable. Most importantly, positive but not 
completely positive maps acting on entangled states may fail to preserve positivity, mapping states 
to non-positive operators. Since such maps will always preserve positivity for separable quantum 
states, which are by definition only classically correlated, positive but not completely positive 
maps can be used to detect the presence of entanglement. The most famous such test is partial 
transposition; the states with no entanglement detectable this way are known as the PPT states 
[25]. In this article, we use free probability and random matrix techniques to construct new families 
of positive but not completely positive maps and apply these maps to the study of entanglement. 

As an illustration of the power of the method, we consider quantum states of a bipartite space 
of finite dimension. Motivated originally by debates about whether room temperature nuclear 
magnetic resonance experiments were capable of producing entanglement, researchers determined 
some time ago that there is a neighborhood of separable states in the vicinity of the maximally 
mixed state [9] , and have produced successively sharper estimates on the size of that neighborhood 
[15]. We consider the related question of perturbing the maximally mixed state by adding to it a 
random Hermitian traceless matrix. This question has been studied in turn by Aubrun, Szarek, 
Ye, and others [2, 3, 4]. In the limit where at least one of the systems has large dimension, one 
finds sharp transitions between separable states, PPT states and, hnally, states whose entanglement 
can be detected using partial transposition. We construct new positive maps capable of detecting 
entanglement to which the partial transposition test is insensitive. 

Formally, let A, B be two C^-algebras. Positive elements in a C* algebra are elements x that are 
self-adjoint and that can be written as x = yy* for some y in the C*-algebra. A positive linear map 
$ : .A —>• 13 is a map that sends a positive element to a positive element. It is well-known that if ^ 
is positive, then = id*, 0 : A4fc(C) 0 ^ A4fc(C) (8) 13 is not necessarily positive, unless A or 

B are commutative. For example A = B = A4n(C),n > 2 and <I> is the transpose map, then, <I>fc 
is not positive as soon as k > 2. 

However if <I>fc is positive, then for I < k, the map is clearly positive too. A map such that 
<I>fc is positive is said to be ^-positive. A completely positive map is one that has this property 
for all integers k. Note that if A. = A4n(C) and B = A4(i(C), min(n, dj-positive is equivalent to 
completely positive [29, Corollary 4.19]. There are many examples of maps that are k positive but 
not {k + l)-positive, but the classification of positive, or even fc-positive maps is far from complete. 

One of the aims of this paper is to provide a new systematic method to obtain fe-positive maps 
from A4ri(C) —)• A4rf(C), based on free probability techniques. More precisely, for all integers n, k 
such that n > k, we describe infinite families of maps $ : Ain(C) —>• A4d{C) that are fe-positive but 
not {k + l)-positive. Our families are defined for any d large enough (and the threshold depends on 
n, k), and their positivity properties depend on the behaviour of the support of the free convolution 
powers of a given compactly supported probability measure. 

Our first main result is Theorem 3.2, and can be stated as follows (we refer to Section 3 for 
details): 

Theorem 1.1. In a non-commutative probability space where a copy of the nx n matrices A4n(C) 
is free from a self-adjoint element X of distribution p,, consider the map : M.n{C) —)• A4 given 
by = EiiXEji, with E = the basis of matrix units. Then, the map <I>^ is 

k-positive if and only i/supp(^®"'/^) C [0, oo). 

The second main series of results looks into a random matrix approximation of the above theorem, 
and focuses on the Gaussian setup. We summarize the results below, and refer to Section 5 for 
details. 
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Theorem 1.2. Let be a shifted GUE random matrix in A4n(C)<SiAid(C) of mean 2 and variance 
a £ [0,1). Then there is a linear map ■ -A4n(C) —)• Aln(C) such that its Choi matrix is again a 
shifted GUE matrix (which depends on Zd), of mean (2 + e)ly/n and variance 1. With probability 
one as d —)• oo; 

(1) The matrix Zd is positive and PPT. 

(2) The map ^d is positive. However, ^d is not completely positive, provided that 2 + e < y/n. 

(3) As soon as 2(2 + e) < a^/n, the map detects the entanglement present in Zd- 

Our paper is organized as follows. Section 2 sets up some notation about positive and fc-positive 
maps, and recalls some results about their Choi matrices. Section 3 is about free probability: it 
recalls abstract facts about freeness, the free additive convolution semigroup, and it gives our main 
new construction of ^-positive maps. Section 4 is about the finite dimensional random approxima¬ 
tion of the model of Section 3, and Section 5 specifies to the case of non-centered GUE random 
matrices, with applications to typical radii of convergence for PPT vs. separability, along with 
non-decomposability results for our random maps. 

Acknowledgements. B.C.’s research was partly supported by NSERC, ERA, and Kakenhi funding. 
P.H.’s research is supported by CIFAR, FQXi and the Simons Foundation. I.N.’s research has been 
supported by a von Humboldt fellowship and by the ANR projects OSQPI 2011 BSOl 008 01 and 
RMTQIT ANR-12-IS01-0001-01. Both B.C. and I.N. were supported by the ANR project STOQ 
ANR-14-CE25-0003. 


2. Positive maps: notations and facts 


We denote by A4n(C) the algebra of ra x re complex matrices. Let E = (Ejj)j , 1 } be the basis 
of matrix units, i.e. EijE^i = Endjk and E*j = Eji. For a C*-algebra A, let : A4„(C) —)• .A be a 
linear map. Its Choi matrix in the basis E, denoted by G# is an element of Ain{C) ® A defined as 

n 

G(j) = ^ ^ Eij (8) 

*j'=i 

Denoting by Bd £ 'Si C'^ the normalized Bell state 


( 1 ) 


Bd 


the Choi matrix can also be written as 






C^ = d[iddSmBdB*d). 

A map $ : M.n{C) —)• A. is called k-positive if the dilated application 

<!> (g) idfe : Mn{C) S Mk{C) —A. (g) Mk(C) 

is positive. Moreover, such a map is called completely positive if it is ^-positive for all k. According 
to a celebrated result of Choi [13], is completely positive if and only if is positive in A4n(C)(gA.. 
A characterization of /c-positivity using the Choi-Jamiolkowski isomorphism has been obtained in 
[26]. 

We recall the following result from [22], characterizing fc-positivity of maps. 

Proposition 2.1. [22, Proposition 2.2] Consider a linear map $ : M.n{C) — )• BiH), where H is a 
Hilbert space. The following are equivalent: 

(1) The map $ is k-positive. 

(2) For all vectors x = Vi S Zi, with yiS Zi £ <CA S H, {x, C^x) > 0. 
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(3) The operator {In®P)C<^{In®P) is positive semidefinite for any rank k orthogonal projection 

PeBin). 

(4) The operator {P (g) (g) 1-^) is positive semidefinite for any rank k orthogonal pro¬ 

jection P G A4„(C) (1-^ is the identity operator in 

3. A MODEL VIA FREE PROBABILITY 

We start this section by recalling some facts from free probability theory that will be used 
throughout the paper. We then introduce our construction of linear maps from probability measures 
and we characterize their block-positivity properties (random matrix versions of these results will 
be considered in Section 4). Finally, we discuss some simple, but illuminating examples in the final 
subsection. 

3.1. Preeness and free products. The monographs [33] and [24] are excellent introductions to 
free probability theory; here, we recall only the concepts needed in this paper. 

A *-non-commutative probability space is a unital *-algebra A endowed with a tracial state f, 
i.e. a linear map 4>: A ^ C satisfying 4>{ab) = (f){ba),4>(aa*) > 0,(()(1) = 1. An element of A is 
called a non-commutative random variable. 

Let Ai,... ,Ak be subalgebras of A having the same unit as A. They are said to be free if for 
all Oj G Aj^ {i = 1,..., k) such that = 0, one has 

4>{ai • • • afc) = 0 

as soon as ji j 2 , j 2 / JSj ■ ■ ■ ,jk-i A Jk- Collections ^i, 52,... of random variables are said to be 
free if the unital subalgebras they generate are free. 

Let («!,... ,afc) be a fc-tuple of self-adjoint random variables and let C(Ai,..., A^) be the free 
♦-algebra of non-commutative polynomials on C generated by the k self-adjoint indeterminates 

Ai,...,Afc. 

The joint distribution of the family is the linear form 

T{ai,...,ak) ■ C(Ai, . . . ,Xk) C 

P 4>{P{ai ,... ,Ofc)). 

Given a /c-tuple (oi,... ,afc) of free random variables such that the distribution of a* is pai-, the 
joint distribution P(ai,...,ak) is uniquely determined by the Pafs. 

In particular, Pai+a 2 Paia 2 depend only on pai and pa 2 - The notations Pai+a 2 = hai S Pa 2 

and Paia 2 = hai ^ Pa 2 Were introduced in Voiculescu’s works [31, 32]; the operations ffl and Kl are 
called the free additive, respectively free multiplicative convolution. Moreover, for any compactly 
supported probability measure p, one can define its free additive convolution powers for any 
real parameter T > 1, see [24, Corollary 14.13]. 

The following result will be quite useful for us and was proved by Nica and Speicher [24, Exercise 
14.21]. 

Lemma 3.1. Let a,p be free elements in a non-commutative probability space {A,(j)) and assume 
that p is a self-adjoint projection of rank t G (0,1) and that a is a self-adjoint random variable 
having distribution p. Then, the distribution oft~^pap in {pAp, 4>{p-)) is p^^P. 

Finally, we discuss free products. Given {Ai,(ji)i^i a family of non-commutative probability 
spaces, there exists a non-commutative probability space {A,4>) in which all {Ai,(f)i) can be em¬ 
bedded in a trace-preserving way such that they are free with respect to each other. This is called 
Voiculescu’s free product construction, and we will need it below. There exists a *-algebra version 
of this construction, a C^-algebra and a von Neumann algebra construction. If A is generated by 
the Aj’s, we write (A, f) = *iei{Ai, (pi). 
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3.2. Maps associated to probability measures. We now introduce the main idea of this paper, 
a construction of a linear map depending on a probability measure We start with a compactly 
supported probability measure // on M and fix an integer n. We will be interested in the space 

/i). Note that this space is spanned as a von Neumann algebra by the operator x ^ x, 
which we will denote by X. By construction, X is a self-adjoint operator and its spectrum is the 
support of n- 

We consider the von Neumann probability space free product [33] 

(7W,tr*]E) := (7Wn(C),tr)*(L°°(M,^i),E) 

and the contracted von Neumann probability space (M, r) where Ai = EnAtEn is the contracted 
algebra of Ai and r is the restriction of tr * E appropriately normalized (by a factor n). 

We consider the map : AiniC) At given by 

^fj.{Eij) = EiiXEji. 

Our main result is as follows: 

Theorem 3.2. Let ^ be a compactly supported probability measure. The map defined above is 
k-positive if and only i/supp(^®"'/^) C [0, oo). 

Proof. The map was constructed in such a way that its Choi matrix has distribution //. Indeed, 
the von Neumann algebra M is a corner of M, so the map a e-)• Eij 0 EuaEji implements 

an isomorphism between Ai and AtniAt); hence, = X. For any rank k orthogonal projection 
P G A4n{C), the distribution of the element is [{l—k/n)6o+k/n5i]Mp,. By the 

Nica-Speicher result 3.1 above, this distribution has positive support if and only if supp(/i®"'/^) C 
[0, oo). The result follows now from an application of Proposition 2.1. □ 

The remainder of the paper is devoted to the study of the maps <1>^ and to applications of the 
above construction to Quantum Information Theory. 

3.3. Examples. Let us start the subsection devoted to examples by making some general purpose 

remarks. First, note that in order for the condition in Theorem 3.2 to hold for some value of k, 
the measure // must have non-negative average. Indeed, the average of is n/fc times the 

average of /r. It follows that, if one starts with a probability distribution p, having negative mean, 
the measures cannot have support in [0,oo). In particular, in order to obtain non-trivial 

examples, we have to rule out any distributions which are symmetric with respect to the origin. 

Similarly, a measure p whose support is already in M+ will not be interesting for us because 
the corresponding map <I>^ will automatically be completely positive and thus /c-positive for all 
\ < k < n (recall that we are ultimately interested in positive maps which are not completely 
positive). In the following, we are interested in non-trivial examples of measures p, i.e. measures 
such that 

supp(/i) ^ [0, oo) but supp(/r®^) C [0, oo) for some T > 1. 

3.3.1. The semicircular case. We start with semicircular distributions, which will play a central 
role in Section 5. As mentioned above, in order to produce non-trivial examples of maps 
one should rule out symmetric distributions p] we consider next shifted (non-centered) semicircle 
distributions. Let us consider SCao- the Wigner semicircle distribution of mean a and variance <7^: 


SC, . = 


V^4(j2 - (x-a)2_^ 

l[a—2cr,a+2cr] 


Note that the support of the semicircle measure is [a — 2(T, a -|- 2 ct]. The semicircle distribution 
plays a role in free probability theory analogous to the role played by the Gaussian distribntion in 
classical probability. Its free cumulants are all zero, except for the first two, which read ki = a, 
K 2 = (t"^. Since the free cnmulants linearize the free additive convolution [24, Proposition 12.3], 
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their values for SCa,"^^ are ki = an/k, respectively K 2 = a’^njk, all the others being zero. Hence, 
SCa,CT^^ = SC ,, /—rr, and thus the support of is 

CL7X j K/ if C7 a / 7% / K 


supp 



an 


- 2a 


n an 

I’T 


+ 2a 


In particular, we have: 

Proposition 3.3. Let n be an integer, a £ M, and a > 0, and consider the map : A^„(C) —)■ 

M associated to the semicircle distribution SCa,o-. The map is k-positive if and only if 

k < a^n/(4(T^). In particular, for any n and any 1 < k < n, there exist parameters a,a > 0 such 
that the above map is k-positive but not k + 1-positive. 


3.3.2. The Marchenko-Pastur case. Here we consider ma,t to be the distribution of the random 
variable 1 — aXt where Xt has a Marcenko-Pastur (or free Poisson of parameter t) distribution 
and a > 0. We recall that the free Poisson law of parameter t is given by (see, e.g. [24, Dehnition 
12 . 12 ]): 


TTt = max(l — t, 0)(lo + 


^4t- (x - 1 - 

2tj-x ^[l+t-2Vi,l+t+2yt]“®> 


and that it has the semigroup property vr^ ffl Tr^ = (this follows from the fact that the free 
cumulants of tt ^ read Kp = t,\/p > 1). In particular, this implies for our purposes that the bottom 
of the spectrum of ma^t is located at 1 — a(l +1 + 2y/i). It follows from Nica and Speicher’s Lemma 
3.1 that if p is a free projection of rank k/n, the bottom of the spectrum of the operator p{l — aXt)p 
is located at 


1 - 


n 


a (^1 + t— + 2y^ tn/k'^ 


We summarize what we obtained above as follows. 


Proposition 3.4. Let n be an integer, a,t > 0, : A4„(C) —A4 whose Choi map is associated 

to the non-commutative random variable 1 — aXt. Ifl — ^a{l-\-t^-\-2-\Jtn/k) > 0 butl —^^^a{l-\- 

+ 2y^ tn/{k + 1)) < 0 then the map is k-positive but not {k + l)-positive. Moreover, for 
each n, k there exists a choice of a, t that satisfies the above inequalities. 


3.3.3. The ‘small rank projection’ case. Now we are interested in the case where X = 1 — aP where 
P is a projection of rank e small (i.e. t{P) = e G (0,1)). 

Proposition 3.5. Hs e —)• 0, the threshold in a to obtain k-positivity converges to n/k. 


Proof. Let us consider some fixed e > 0. Then, there exists a variable Y having a Marchenko-Pastur 
distribution of parameter e such that, for any rj > 6y/e, 


1 — (a + r])Y < X < 1 — {a — r])Y. 


Indeed, recall that the Marchenko-Pastur distribution of parameter e < 1 can be written as 


TTe = (1 - s)6o+eTTe, 

where is a probability measure supported on the compact interval [(1 — (1 + \/e)^]. The 

result follows after applying Proposition 3.4 and letting first e —)• 0 and then rj ^ 0. □ 
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3.4. Quantum states detected by “free” maps. Let be the set of n-dimensional quantum 
states 

Vn := {p S Mn : /O > 0 and Tr/? = 1}, 
and consider the set of separable states 

<S£T^n^m — \ ^ ^ tiPi ® . ti ^ 0, Pi ^ T-^ni ^ ^ f 

t i i 

Given a positive map : A4n —)• Md, we introduce the convex body of n • m-dimensional quantum 
states which remain positive under the partial action of 

— {p G T-^nm ■ ® idm] (p) ^ 0}. 

We first record the following elementary result: 

Lemma 3.6. For any positive map $ and any dimension m, the set is a convex body and it 

contains the set S£Vn,m of separable states. If 4* is not m-positive, then K^^rn is a strict subset of 
T> 

^nm • 

Let us turn our attention to the particular case of a map 4*^ as defined in the previous sections. 
We start with the following proposition, showing that any non-trivial detects pure entangled 
states. 

Proposition 3.7. Let p be a compactly supported probability measure such that supp(p®"') C [0, oo), 
but supp(p®”/^)n( —oo, 0) / 0, for all 2 < k < n. Then, the set contains all separable states 

but no pure entangled state. Put differently, <1>^ detects perfectly the entanglement of pure states. 

Proof. The claim about separable states follows from the fact that the map is (l-)positive. Let 
p be a rank one self-adjoint projection on A4„(C) (8) A4fc(C) that is not separable (i.e. p is a multiple 
of an entangled state). Let q G A4n(C) be the orthogonal projection on the range of the partial 
trace of p; since p is entangled, rk(g) > 2. 

By direct inspection, studying the positivity of [<1>^ (g) idm](p) is exactly the same problem as 
studying the positivity of (g® l;v()X(g(8) 1 _a 4); indeed, using a Schur map argument, one can assume 
without loss of generality that the eigenvalues of the partial trace of p are all equal. However, by 
our assumptions, the latter element does not have a positive spectrum. Therefore, does not 

contain pure entangled states. □ 

Let us now introduce two intermediate sets between the set of separable states and the set 
of all quantum states, corresponding to the reduction criterion [10, 20] and the positive partial 
transposition (or PPT) criterion [25]: 

TZ£'Dn,m = {/O G F>nm ■ In ® [Tr^ (g idm](/0) “ P > 0} 

PPT n,m = {/O G Vnm : [T„ (g) id^] {p) > 0}, 

where T denotes the transposition map. It is known [10] that these sets satisfy the following chain 
of inclusions (which all become strict for large enough values of n, m) 

S£Pn,m T PPTn,m P£Pn,m T Pnm- 

We show next that the maps corresponding to the “small rank projection” measures (see Section 
3.3.3) detect the same entangled states as the reduction criterion, in a suitable asymptotic regime. 

Theorem 3.8. In the setup of subsection 3.3.3, with parameters a = n, e —)• 0, the convex body 
K^,m tends to P£Vn,m- 

Before we prove this result, let us provide an equivalent, but more general form of the reduction 
criterion. 
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Lemma 3.9. Let $ : Ain(C) —>■ A be a (non-unital) embedding of *-algebras, sending the matrix 
units Eij of Ain{C) to elements Fij G A. Then, for any integer m and for any quantum state 

n n 

In ® [Tr„ (g) idm](/o) - /9 > 0 (8) ^ Pii - ^ Fij ® pij > 0. 

2 = 1 hj = ^ 


Proof of Theorem 3.8. As in the context of section 3.3.3 we are looking at a Choi matrix of the 
form = 1 — nP where P is a projection of trace e. Let be the corresponding map (obtained 
from via Theorem 3.2). 

For a given T>nm ^ 9 = Y17j=i ® PL ^ '^nm, we have 

n n 

(g) idm](p) = Im <8 ^Pm - ^ Pj;-^ (g) pij > 0, 
i=l *,i=i 

where pj^j^ = nEuPEji. Let us show now that the pj-j'^ are approximate matrix units, i.e. they 
verify the ^-algebra relations for matrix units up to some vanishing precision as e —^ 0. To show 
that, as e —0, for all i,j, k, I, 

( 2 ) 

we just have to compute 

n^EuPEjiEikPPii - nSjkPiiPEii = nEu{nPEjkP - SjkP)Eii. 

We show next that, as e ^ 0, 


(3) \\nPEjkP — djkP\\ —)• 0. 

Let us first consider the case j = k. Again from Nica and Speicher’s Lemma 3.1, it follows that the 
distribution of the element P{nEjj — I)P in the reduced non-commutative probability space PAiP 

is that of the dilation of the free additive power Pg (^n~^Sn-i + (1 — . Let e = 1/A^, 

for some (large) integer N. Then, by the super-convergence result of Bercovici and Voiculescu [7, 
Theorem 7] (see also [23, Example 1]), we have that the support of 

, fflWl 


Dm- 


Af-l/2 


(n ^6n-i + (1 - n ^)5-i) 


is contained, for N large enough, into some fixed compact interval. The normalization in N~^ 
shows the desired convergence in norm to zero. The case j A k is treated in a similar manner, by 
considering the symmetrized version \\P{Ejk -|- Ekj)P\\ = 2\\PEjkP\\. 

By multiplying on the left and on the right equation (3) by Eu, resp. En which are both of 
operator norm one, one obtains (2); note that in all the computations above, n is a fixed parameter. 

It follows from [8, Proposition 11.8.3.22] that approximate matrix units can be approximated by 
true (exact) matrix units: there exist true matrix units G Ai such that, for all i,j, — 


^(e)| 




0 as e —)• 0. Define 


^e,m.{p) = i-M 

By the previous lemma, p G IZ£T). 


n 




Pij P 0. 


2=1 = l 

iff. 'I'£,m(p) > 0 in A4. But we have 


||[^>g (g) idm](p) - ^£,m(p)|| = 
finishing the proof. 




- G^f)) ® p,j 




< sup \\Fj;f - G[f\ 


0 , 


□ 
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In other words, an infinitesimally small free projection captures exactly the set of quantum states 
whose entanglement is detected by the reduction criterion. Since TZEVn^m 3 'PVTn,m-, an important 
question would is to check whether the sets can be compared to WTn,m for appropriate 

measures /r. Although we do not address this question directly, we will see with Theorem 5.4 that 
with random techniques, a result of this flavor holds. 

We finish this section by a general open question. For all n, k, I, we introduce the set 

KlTn, ■= n 

fj,: jc[0,oo) 


We call it the collection of freely detectable quantum states in T>nm- 

In view of the preliminary results of this section, we believe that this set is of some interest in 
quantum information theory. The sequel of this paper will show that suitable approximations of 
this set have a very strong entanglement detection power. This leads us to stating the following 
question, that has to our mind an interest both in free probability and quantum information theory: 

Question 3.10. Can one give a description of We know this set contains the collection 

k-separable states (Theorem 3.2), but how much bigger is it? Are there values of the parameters 
n,m,m for which the set is precisely the set of k-separable states from Vnm? 

4. A FINITE DIMENSIONAL, RANDOM VERSION 

Theorem 3.2 and the construction of the maps give a way to construct new families of positive 
but not completely positive maps. The results in Section 3 are not very concrete in the sense that 
the compression of free products of von Neumann algebras is a rather abstract object. In this 
section, we show that it is possible to interpret the previous results as the limiting case of random 
maps between finite dimensional matrix algebras. In particular, we obtain a large class of linear 
maps between matrix algebras that are k positive but not A: + 1 positive. For this we need a result 
about norm approximations. 

4.1. Norm behaviour of products of random matrices. We recall that if A is a d-dimensional 
self-adjoint matrix, its empirical eigenvalue distribution is d~^ where Aj are the eigenvalues 

of X. For any probability measure p. on the real line, its repartition function is defined as : t 
/i((-oo,t]). 

For the purposes of this section, we will say that a sequence of repartition functions fd converges 
to a repartition function / if for all e > 0, there exists an do such that for all d > do, 

Vt e M, f{t-e) - e < fd{t) < f{t -h e) + e. 

Theorem 4.1. Let Ad,Bd he independent positive self-adjoint random matrices in Aid{C), such 
that at least one of Ad or Bd has a distribution invariant under unitary conjugation. Let fd be 
the repartition function of Ad and gd he the repartition function of Bd. Assume that the (a priori 
random) repartition functions fd,gd converge almost surely respectively to f,g which are repartition 
functions of two self-adjoint, bounded and freely independent random variables x and y. Assume 
also that the operator norm of Ad (resp. Bd) converges to the operator norm of x (resp. y). Then, 
almost surely as d —)■ oo, 

\\AdBd\ \ ||xy||. 

The above result is a special case of a much more general class of results, see [11, 5]. 
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4.2. The finite dimensional, random matrix model. In this section, we fix again a compactly 
supported probability measure ^ and an integer n. For each integer d, we introduce a real-valued 
diagonal matrix G A4n(C) (8i Ald(C) whose empirical eigenvalue distribution converges to /i and 
whose extremal eigenvalues converge to the respective extrema of the support of fi: 

lim Amin(-’^d) = min{x : x G supp(^)} 

d^OQ 

lim Amax(-’^d) = max{x : x G supp(/i)}. 

d—^oo 

Let also Ud be a random unitary matrix in lAnd distributed according to the Haar measure on the 
unitary group. We introduce the map 4)^^^ : -Adn(C) —)• A4d(C) whose Choi matrix is the random 
matrix Cd '■= UdXdU^ G X4n{C) (8) Ald(C); this can be done using the hermiticity preserving map 
GA self-adjoint matrix incarnation of the Choi-Jamiolkowski isomorphism. More precisely, we have 

VX G Mn{C), ^^,d{X) = [Tr ® id](C7rf • ® I^). 

We can now state the main result of this section. 


Theorem 4.2. The sequence of random linear maps {^^^d)d has the following properties: 

(1) //supp(/i®"'/^) C (0, oo), then, almost surely as d ^ oo, is k-positive. 

( 2 ) //supp(/i®’^/^) n (— 00 , 0 ) / 0, then, almost surely as d ^ oo, ^^^d is not k-positive. 


Proof. Let Cd be the Choi matrix of the map ^fj,^d, he. Cd = UdXdU^ for a sequence of real 
diagonal matrices {Xd)d converging in distribution to p,. The assumption implies that there exists 
some positive constant (5 > 0 such that supp(//®”/^) C [(j, oo). Fix M > maxsupp(/x) and let e > 0 
be a small constant which will be fixed later. Fix an e-net Q of projections of rank k on A4n(C) 
(for the operator norm). This net can be chosen to be finite because the collection of projections 
is a compact space. In other words, for any rank k projection P G A4n(C), there exists a rank k 
projection Q G Q such that \\P — Q\\ < £. 

For a fixed element Q G Q, we obtain, using Theorem 4.1 and Lemma 3.1, almost surely as 
d —)• oo. 


lim An 
d—)-oo 


{Q®Id-Cd-Q® supp(/x®’"/^) > 


k 


kS 


n 


Taking the intersection of a finite number of almost sure events, we obtain the same conclusion for 
all elements Q oi Q. Moreover, since, for all rank k projections P, we have almost surely 


P^Id-Cd-P®Id-Q^Id-Cd-Q^Id\\<\\P-Q\\- \\Cd\\ < Ms, 


we conclude that, almost surely as d —)• oo. 


lim min Amin 

d—^oo P projection of rank k 


{P ® Id - Cd - P ® 


>- Ms, 

n 


quantity which is positive as soon as e < {k5)/{nM). By Proposition 2.1, this completes the first 
part. 

To prove the second point, consider any fixed orthogonal projection P G A4n(C) of rank k. Using 
Theorem 4.1, we obtain that, almost surely as d —)• oo, the smallest non-trivial eigenvalue of the 
matrix {P (8) Id)Cd{P <8> Id) converges to k/n times the bottom of the support of which is 

a negative constant. Therefore, by Proposition 2.1, the map is not positive, almost surely as 
d —)• oo. □ 


When comparing the result above with Theorem 3.2, the reader will notice that one needs a 
gap around zero in the random matrix formulation. This is due to the fact that one needs to 
avoid eigenvalues fluctuating around zero in order to draw a conclusion about the positivity of 
a matrix. Finally, from the above theorem, one can immediately deduce the following corollary, 
simply by replacing the almost sure notion of convergence of random variables with the weaker one 
of convergence in probability. 
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Corollary 4.3. The sequence of random linear maps has the following properties: 

(1) // supp(^®"'/^) C (0, cxd), then 

lim P is k-positive] = 1. 

d^oo ’ 


(2) //supp(/i®"/^) n (— 00 , 0 ) / 0, then 

lim P 'Is k-positive] = 0. 

d^oo ’ 


5. GUE QUANTUM STATES AND MAPS 

In this section, we focus on (possibly traceless) random matrices from the Gaussian Unitary 
Ensemble (GUE). The second chapter of the monograph [1] is an excellent exposition of the subject. 
The interest of GUE random matrices is two-fold. Firstly, it is the case where the norm behavior 
of random matrices is best understood. Secondly, GUE matrices are both unitarily invariant (their 
distribution is invariant with respect to unitary conjugations) and Wigner matrices (the entries 
of the random matrix are independent and identically distributed, up to symmetries). We write 
i = ■\/—l in order to avoid confusions with the index i of matrix elements. 


Definition 5.1. A random matrix Z E Alrf(C) is said to have the GUEj^ distribution if its entries 
are as follows: 


Zij — < 


^iij ^fd-j 
{Xij iYij) 


Z 




if i = j 
if i < j 
ifi > j, 


where {Xij,Yij}fare i.i.d. centered, standard real Gaussian random variables. 

We recall next the Wigner theorem, as well as the stronger, norm convergence of GUE matrices. 


Theorem 5.2. Consider a sequence of random matrices X^ G Ald(C). Then, almost surely 

as d ^ oo, the empirical distribution of X^ converges weakly to a standard semicircular distribution: 

1 

- J]Ai(Xrf)^SCo,i. 

i=l 

Moreover, the norm of X^ converges almost surely to the maximum of the support: 

lim IIX^II = 2. 
d^co 

We start with the following simple but important lemma, which will allow us to construct different 
realizations of GUE matrices acting on a tensor product space. Below, that the matrices {Eij} 
denote the matrix units of Alri(C). 


Lemma 5.3. For any integers n, d, consider a family of i.i.d. GUE^ random matrices {Xij 
Define 


Mnd(C) BZ:= 



n 

'y ^ Eii 0 Xii 

i=l 


+ 



'y ^ Eij <8* {Xij Y iF)j) 


+ 



^ ^ Eji < 8 ) {Xij iYijf 

l<i<j<n 


Then, Z has a GUE^^ distribution. 
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Proof. The (nd)^ elements of Z are centered, Gaussian, independent (up to symmetries) random 
variables. The normalizations are chosen in such a way that the variances are exactly those of a 
GUE„rf random matrix. □ 


5.1. PPT entanglement. We are now going to construct an asymptotically positive semidefinite 
and PPT matrix G Aln(C) (8) Alrf(C) as well as a map which will detect the entanglement in 
Zd- Both constructions are inspired from Lemma 5.3. 

Gonsider a fixed integer re > 1 and fixed parameters e > 0 and a G (0, 1). Let and 

V‘^ij}i<i<j<n be two families of i.i.d. GUE^ random matrices. Define the matrices 


(4) 


and 


Mnd{C) 3 Zd := 2In ® Id — a 



n 

^ ^ Eii ® Xii 

i=l 


+ 



^ ^ Eij ® {Xij + ) 


1 



^ ^ Eji ® {Xij 

l<2<ji<n 


2 1 ^ 

Xind{II) 3 Cd •= —® Id~\ - 1= ^ ^ Ejj ® W 

■ore -ore —' 

1=1 


(5) 


+ 


+ 


’ \ < 't ^ rt < 'rt 


l<2<j'<n 






l<2<j'<n 


and consider the linear map : Ad„(C) —)• Aid{C) whose Choi matrix is Cd- In other words, for 
an input p G M„(C), 

n 

( 6 ) = Y 

*J=1 

where Cd{i,j) G Add(C) denotes the {i,j) block of Cd- 


Theorem 5.4. For the linear maps defined in (4),(5), and (6), the following statements hold: 

(1) For all values ofn, a, almost surely as d ^ oo, the random matrix Zd is positive semidefinite 
and PPT. 

(2) For all values of n,e, almost surely as d ^ oo, the linear map is positive. 

(3) If Bd denotes the normalized Bell state (1), then, almost surely as d ^ oo, 

lim {Bd, [4>d ® idd]{Zd)Bd) = 2(2 + £)y/n - an. 

d—^oo 

In particular, for any choice of parameters such that 2(2 + e) < a^/n, the matrix Zd is, 
almost surely as d ^ oo, asymptotically PPT and entangled. 

Proof. Eor the first point, note that the matrix Zd can be written as Zd = 2Ind — ctZd, where Zd is 
a GUE^rf matrix (see Lemma 5.3. It follows that the support of its asymptotic spectral distribution 
is the interval [2 — 2a, 2 + 2a], which is supported on the positive half-line; strong convergence of 
Zd proves that Zd is asymptotically positive semidefinite. It is also PPT, since transposed GUE 
matrices are still GUE. 
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For the second claim, to show positivity of the map having as a Choi matrix, we are going 
to use Theorem 4.2. Indeed, at fixed n, as d —)• oo, the random matrix Cd has asymptotic spectral 
distribution jj = SC( 2 +e)/y 7 j^i. We get /x®”' = SC( 2 +£).yH,v^) which has support [s^/n, (4 + e)\/n] C 
(0,oo), hence ^d is asymptotically positive, by Theorem 4.2. 

Let us now show that the third point holds. Denoting by A{i,j) G A4d(C) the (i,j) block of a 
matrix A G A4„d(C), a direct linear algebra computation shows that 


{Bd,[^d®^dd]{Zd)Bd) = -,Y. Tr[Zrf(z,j)Q(z,j)®] 




2=1 




Using, for all i,j, almost surely as d —>• oo 


and 


lim -Tr(Xf ) = lim -Tr{Y^A = 1 
d^oo d ^ d^oo d ^ 


lim ^YviXijYij) = lim ^TA{Xij) = lim \Tv{Yij) = 0, 
d^oo d d-^OQ d d^oo d 


we can conclude. 


□ 


The construction of the matrix Zd can be interpreted as follows. Let us replace the notation Zd by 
Zd{a), i.e. view the Gaussian component of Zd as a ‘direction’ and a as an intensity when choosing 
at random from the vicinity of the maximally mixed state (everything up to trivial scalars). With 
probability 1 we have that as d —)• oo, Zd{B) is a multiple of a PPT state, and that actually it is 
on the boundary of PPT in the sense that if a = 1 is replaced by a greater value of a, it will leave 
the convex cone of PPT matrices. 

The following corollary is itself of interest: 

Corollary 5.5. For any n > 16, there exist a,£ such that the map detects the entanglement in 
the PPT matrix Zd- 

In other words, as soon as n > 16, the maps we introduced in sections 3 and 4 can detect PPT 
states that are not separable. The analysis above also shows that the maps become very efficient 
as entanglement detection tools for n large. Indeed, 

Corollary 5.6. For any 1 > a> y/n the positive matrix Zd is PPT entangled. 

This following result is a direct extension of Theorem 5.4. 

Proposition 5.7. Fix parameters n,l such that n > 16/. For any 1 > a > Ay/lly/n, the positive 
matrix Zd from (4) is PPT but not l-separable, with high probability as d ^ oo. 

Proof. The proof goes exactly as for Theorem 5.4, except that one has to replace positivity by 
/-positivity, which affects the threshold accordingly. □ 


5.2. Random separability. Finally, we would like to finish this section by mentioning that the 
threshold 4/-y/n from Corollary 5.6, although perhaps not optimal, is of the right order with respect 
to n, when n is large. 
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Theorem 5.8. Let he a standard GUE random matrix in Aindi^)> a G (0,1) he a fixed 
constant. With probahility 1 as d ^ oo, the matrix 2Ind+ctXd is separable as soon as ^ . 

In particular, if a < y/nj{2(n — 1) + y/n), the hermitian matrix defined in (4) is, almost surely 
as d ^ oo, separable. 

Proof. Consider a fanily {xj;j\ X^J\ X^p, xj;j^}i<i<j<n of i.i.d. GUE^ operators. Eor /3 G (0,1), 
define, for 1 < z < j < n, the operators 

= (En + E,, + Ej, + E,j} ® ifixjp + 21,) 

= {Eu - E^j - Ej, + E,fi ® {fix\f + 21,) 

T^f = {Eu + iE,j - iEj, + E,j) ® + 21,) 

T^f = {Eu - iEij + iEj, + Ejj) ® (pX^f + 21,). 

Eor 1 < i < n we further introduce 


Ti = Ei, 


E (Ei ^ + E? + Ef + Ef) + 2/3^'^ - 8(n - 1)1, 


j¥=i 


where we write, for convenience, X^-^ = X!>I for i < j. Let 


X) 






Y,= 


2l3y/n 


E ELf+Et 

l<i<j<n s=l 


i=l 


A direct computations shows that Y, follows a GUE^d distribution. Define 

4 n 


(7) 


Y, . xl^i, 4“ Y, 


1 


2f3y/n 


E Y'dI + Y.r 

l<i<j<n s=l 2=1 


where the new diagonal blocks read 


Ti = Ei, 


E (Ei ^ + Ef + Ef + Ef) + 2/3E'^ - 8(n - 1)1, + 2xfiV^I, 


j¥=i 


We look now for the range of x such that equation (7) is a separable decomposition for the matrix Y,. 
Since /3 < 1, each of the operators fixjp + 21, is, almost surely as d ^ oo, positive semidefinite; 

hence, the same hols for all Moreover, for all i, the limiting eigenvalue distribution of the 

random matrix T is 


SC 

which has positive support if and only if 


2x^y/n—S{n—l),2l3y/n^ 


x>2 + 


4(n — 1) 


fiVn 

Optimizing over the values of /3 G (0,1), we conclude that xind + GUE„rf is separable whenever 

4(n — 1) 

X > 2 + -^ ^ 


n 


Eor the matrix Z, in (4), we have now that 

Z, = a(^^ + GUE^rf^ 
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is separable whenever 


which finishes the proof. 


2 


- > 2 + 
a 


4(n — 1) 

y/n 


□ 


Corollary 5.9. In almost all directions, for any n > 16, the convex body of separable states is 
contained between two euclidean balls of relative radii 4:{y/n + 2(n — l))/n. For n large, the ratio 
that we can obtain is of the order of 8. 

Note that the comparison is quite interesting with [15, Corollary 3] that says that the largest 
euclidean ball inside separable sets in our above setup is of radius 1/n. Obviously our order of 
^l\fn is much larger, however, our result is true ’only’ in ’almost all’ directions, so it is actually 
compatible - and of intrinsic interest. 

Our results also complement the following finding of Aubrun, Szarek and Ye in [3], which can be 
stated informally as follows: 


Theorem 5.10. Let p £ AinmiC) be a random quantum state distributed along the induced measure 
with parameters {nm,s). Then, there exist constants ci ^2 cmd a function sq = so{n,m) such that 

cinm min(n, m) < sq < C 2 nmlog(nm)^ min(n, m) 

and 

(1) If s < {1 — e)so, the probability that p is separable is exponentially low. 

(2) If s > {1 + e)so, the probability that p is entangled is exponentially low. 


Remark 5.11. In the vein of the above theorem, it would be interesting to find good examples of 
typically l-separable non-centered GUE. Our construction above is specific to the one-separable case, 
and we do not see how to extend it in a non-trivial way to the l-separable setup. 


5.3. Indecomposability. We recall [29, Definition 1.2.8] that a positive map 4> is called decom¬ 
posable if it can be written as 

4> = $1 + T o 4>2, 

where T is the transposition map and $ 1^2 are completely positive maps. 

Theorem 5.12. In the setting of Theorem 5.4, as soon as 2(2+e) < a^/n, almost surely as d ^ 00 , 
the positive linear map having Choi matrix C^, is indecomposable. 


Proof. This follows from the third point in Theorem 5.4: the map detects (for d large) the 
entanglement in the PPT state Z^, hence it cannot be decomposable, see [29, Corollary 7.2.5 and 
Proposition 7.3.6]. □ 
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